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\^ . We have employed time-dependent local-spin density theory to analyze 

^^ \ the multipole spin and charge density excitations recently found in GaAs- 

AlGaAs quatum dots [C. Schiiller et al, Phys. Rev. Lett 80, 2673 (1998)]. 
^ ' The overall agreement between theory and experiment is good, identifying the 

(y-v ' angular momentum of the modes observed in the experiment. We have found 

C^ . that high multipolarity spin density edge modes originate from interband 

transitions instead that from intraband transitions, as it happens in the dipole 



0\ ' case. 
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The characteristic single particle and collective excitations of typical quantum dots (QD) 
are known to lie in the far-infrared (FIR) energy region, i.e., they have energies that, de- 
pending on the size of the dot, span the range from a few tens of meV to a fraction of meV. 
Experimental information about FIR spectra was first obtained from photon absorption ex- 
periments on InSb and on GaAs quantum dots [|l],H. Since the confining potential for small 
dots is parabolic to a good approximation, and in the FIR regime the dipole approximation 
works well, the absorption spectrum is rather insensitive to the number of electrons in the 
dot, measuring to a large extent only the center-of-mass excitations, which at non zero mag- 
netic fields (B) correspond to the two allowed dipole transitions arising from each of the two 
possible circular polarizations of the absorbed light. Two limitations of the absorption pro- 
cess, namely that it is dominated by the L = 1 multipole of the incoming electromagnetic 
wave, and its insensitivity to the electronic spin degree of freedom, have motivated that 
theorists have been mostly concerned with the study of dipole charge density excitations 
(CDE), although calculations of higher multipolarity density modes exist in the literature 
(see for example Refs. [§]-|§]). 

The situation is changing with the use of inelastic light scattering to experimentally 
study QD excitations l^-^ . In a way, these studies complement the similar ones carried out 
in the past on the two dimensional electron gas (2dEG) ||10|- Besides opening the possibility 
to study the wave vector dispersion, inelastic light scattering allows to disentangle charge 
from spin density (SDE) and single particle excitations (SPE), and to observe them all in 
the same sample. 

Very recently, the B dispersion of CDE's and SDE's of different multipolarities has 
been experimentally determined in GaAs-AlGaAs quantum dots |jll|. We present here a 
theoretical interpretation of these results based on the time-dependent local-spin density 
theory (TDLSDT) which addresses for the first time the description of high multipolarity 
spin modes. 

To this end, we have obtained the ground state (gs) of an A^ = 200 electron dot confined 
by a uniform, positively charged disk oi R = 120 nm solving the appropiate Kohn-Sham 
(KS) equations. These values correspond to a quantum dot thoroughly studied by Schiiller 
et al []rT|. The exchange-correlation energy density SxciPji^), where p is the electron den- 
sity and m the spin magnetization, has been constructed from the results of Ref. |jl2| on 
the nonpolarized and fully polarized 2dEG using the two dimensional von Barth and Hedin 
prescription |13| to interpolate between both regimes. The only free parameter in the cal- 
culation is the number of positive charges in the disk, which has been set to A^"^ = 404 
to reproduce the dipole SDE at B = 0. The range of B values investigated in this work 
corresponds to filling factors larger than 3. 

Once the KS gs has been worked out, we have determined the induced densities originated 
by an external multipole field employing linear-response theory. Since we have described 
at length the dipole longitudinal response in dots [Q, we give here only a few details for 
presentation purposes. For independent electrons in the KS mean field, the variation Sp^^^ 
induced in the spin density p^ (a =t, i) by an external spin- dependent field F, whose 
non-temporal dependence we denote as F = J2a fa{r) |o")(cr|, can be written as 



where Xal' i^ ^^^ -^^ s\)m. density correlation function for independent electrons. In this 
limit, the frequency uj corresponds to the harmonic time dependence of the external field F 
and of the induced 5p^^\ Eq. (|l|) is a 2x2 matrix equation in the two-component Pauli space. 
In longitudinal response theory, F is diagonal in this space, and its diagonal components 

are written as a vector -F = <• • We consider the external L-pole fields 

which cause, respectively, the charge and spin density L-modes. For the monopole L = 
mode, these fields are simply taken proportional to r^. To differentiate the induced densities 
of each excitation channel they will be labelled with an additional superscript as 5p^'^^ or 

The TDLSDT induced densities are obtained solving the equation 

Spl''\f,u;) = Sp^^'^\f,uj) + Y: Jdndf,x'^}^if,n;u;)K,,^M,r2)Sp^^\f,u) , (3) 
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where either A = p oi A = m, and the kernel Kcra' {r, r ') is the residual two-body interaction. 
Equations (j^) have been solved as a generalized matrix equation in coordinate space. 
Taking into account angular decompositions of Xaa' and K„„/ of the kind Kf^^/ (r , r ') = 
J2i Kali{r, r')e^^'^^~^'\ it is enough to solve this equation for each multipole separately because 
only modes with i = ±L couple to the external L-pole field. One has 
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where Kn{x) is given by the hypergeometric function |T^ ^F(l/2, n + 1/2; n -|- 1; x^). 

For a polarized system having a non zero magnetization in the gs, the ±L modes are 
not degenerate and give rise to two excitation branches with AL^ = ±L, where L^ is the gs 
orbital angular momentum. The induced charge or magnetization densities corresponding 
to density and spin responses are given by ^p^"^^ = 5p\ + 5p\ and 5rrS^'> = 5p\ — 5p\ . 
From them, the dynamical polarizabilities in the density and spin channels are respectively 
given by 

app{^,u)= I drr\^\+^5p^f\r) 
ammii,uj) = Jdrr\'\+'5m^^\r). (5) 

For each L value, taking into account both ±L possibilities we define a)^AA — c^AAiL, uj) + 
aAAi~L,uj). Their imaginary parts are proportional to the strength functions SAiA = 



Figures m and ^ represent the spin and charge strength functions for L = and 2, respec- 
tively. In the quadrupole case we have indicated with a +(— ) sign the SDE's arising from 
the +L(—L) component of the F^"^^ operator in Eq. ^ (the low energy CDE is always a + 
type excitation, whereas the high energy CDE are — type excitations arising from the cor- 
responding component of F^^'). We have found that the spin peaks are rather fragmented, 
especially in the monopole case. However, they still are collective modes, with energies red- 
shifted from the free electron ones due to the attractive character of the exchange-correlation 
vertex corrections. 

We would like to draw the attention to the — type, low energy quadrupole SDE which is 
seen in Fig. Q to carry an appreciable strength at cu ~ 3.1 meV and B = 2 T. When a mag- 
netic field is perpendicularly applied to a QD, it is well known that the low energy excitation 
modes in the density channel are dipole edge CDE's arising from intraband transitions while 
bulk, interband transitions lie higher in energy. That may change with increasing L, and it 
is particularly easy to see that this is the case if one looks at the SDE's. An inspection to the 
KS single electron energies shown in Fig. ^ reveals that at high L's, interband electron-hole 
excitations are at lower energies than intraband ones. Since the residual electron-hole inter- 
action is weak in this channel, we have found that at B = 2 T, the lowest energy octupole 
SDE is indeed a mode built from interband electron-hole excitations. Still, this is an edge 
mode, as its existence is only possible because of the finite size of the system. 

Figures ^ and ^ display the B dispersion of the more intense SDE's and CDE's, respec- 
tively. The solid symbols represent the experimental data [0]. It can be seen from these 



figures that the overall agreement between theory and experiment is good. In both spin and 
charge density channels, TDLSDT reproduces the weak B dependence of the L = mode 
found in the experiment at small B values. Our calculations confirm the L = 0, 1, and 2 
multipolarity assigned in the experiment to the lower SDE's, but cannot identify the origin 
of the higher SDE. Ruling out the possibility that it is an L = 3 or 4 SDE (see Fig ^), it 
might correspond to a L < 3 mode which gets some strength when the dot is probed by an 
excitation operator carrying a finite momentum q on the dot plane, such as a Bessel function 
JiiQr) instead of the r^ multipole we have been using. Indeed, the signal of that peak is 



weak and broad, as mentioned in Ref. ||TT |. 



At B = 0, the energies of the L > spin density excitations follow the simple rule 
El ~ LEi . We attribute this to the weakness of the residual interaction in the spin channel. 
The prominent role played by the strong residual interaction in the charge density channel 
causes that rule to fail for CDE's. 

As a general trend, the strength carried by the positive B dispersion branch correspond- 
ing to the high L spin density excitations diminishes as B increases. We have also found 
that the spin strength becomes more fragmented with increasing L, whereas the bulk and 
edge magnetoplasmons associated with the ±L excitations are well defined modes. 

The positive B dispersion branches of the CDE's reveal a complicated pattern at inter- 
mediate B values, quite different from the expected classical one holding up to i? ~ 2-3 T, 
but that however fits a large set of the experimental modes. The behavior of these branches 
has an interesting quantal origin, namely the formation of well defined Landau bands for 
magnetic fields larger than a critical value. Above it, the more intense high energy collective 
peaks mostly arise from transitions between Landau bands whose index M differs in one 
unit, AM = 1. Since these bands are made of many single electron states with different i 



values and energies rather £ independent if B is high enough |T6[ , this explains the otherwise 



striking quasi L-degeneracy of the plasnion energies, only broken by finite size effects and 
the L dependence of the residual interaction. Other modes with AM = 2 build branches 
satellite of those formed by the more intense L-peaks, which are clearly seen in the calcu- 
lation. Satellite branches of this kind appear even in the dipole case PJT7[], and are a clear 
signature of nonparabolic confinement. 

This work has been performed under grants PB95-1249 and PB95-0492 from CICYT, 
Spain, and 1998SGR00011 from Generalitat of Catalunya. A.E. and M. B. (Ref. PR1997- 
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FIGURES 

FIG. 1. Monopole strength function in arbitrary units as a function of energy. The thick sohd 
hne represents the charge density strength, the dashed line the spin density strength, and the thin 
sohd line the free electron strength. 

FIG. 2. Same as Fig. [l| for the quadrupole mode. 

FIG. 3. Single electron energies as a function of orbital angular momentum i for i? = 2 T. 
The horizontal line represents the electron chemical potential. Full, upright triangles correspond 
to o =1 states, and the empty, downright triangles to o =J, states. Interband and intraband 
transitions with A£ = 2, 3 and 4 are represented to illustrate the energy crossing discussed in the 
text. 

FIG. 4. Energy of the more intense SDE's as a function of B. The lines are drawn to guide the 



eye, and the solid symbols represent the experimental data [11 1. 



FIG. 5. Same as Fig. for the more intense CDF's. 
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